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The main objective of this work is minimizing the cost of electricity of solid oxide fuel cell stacks by 
decelerating degradation mechanisms rate in long term operation for stationary power generation appli- 
cations. The degradation mechanisms in solid oxide fuel cells are caused by microstructural changes, 
reactions between lanthanum strontium manganite and electrolyte, poisoning by chromium, carburiza- 
tion on nickel particles, formation of nickel sulfide, nickel coarsening, nickel oxidation, loss of conductiv- 
ity and crack formation in the electrolyte. The rate of degradation mechanisms depends on the cell 
operating conditions (cell voltage and fuel utilization). In this study, the degradation based optimization 
framework is developed which determines optimum operating conditions to achieve a minimum cost of 
electricity. To show the effectiveness of the developed framework, optimization results are compared 
with the case that system operates at its design point. Results illustrate optimum operating conditions 
decrease the cost of electricity by 7.12%. The performed study indicates that degradation based optimiza- 
tion is a beneficial concept for long term performance degradation analysis of energy conversion systems. 
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1. Introduction 


In recent years, the interest in electrochemical power genera- 
tion devices based on solid electrolytes has dramatically grown. 
This concern is due to the advantages of solid electrolytes over liq- 
uid electrolytes. For instance, solids can operate in a wider range 
and therefore, have a lower activation loss |1]. However, the long 
term operation with respect to contaminants and degradation 
mechanisms is still mitigates a broad implementation of SOFCs. 
Thus, a fundamental understanding of degradation mechanisms 
to analysis the long term operation of SOFCs is of great interest | 2]. 


1.1. Concept 


In long term operation, solid oxide fuel cell characteristics and 
performance deteriorates dramatically. Operating conditions affect 
the rate of degradation mechanisms. To study SOFC in long term 
operation, considering degradation mechanisms is unavoidable [3]. 

The degradation based optimization (DBO) is an optimization 
framework that considers system degradation mechanisms in the 
optimization procedure [4]. The flow diagram of the developed 
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DBO framework for SOFC stacks is presented in Fig. 1. At the initial 
time, system characteristics are known. In the degradation model, 
the rate of degradation can be calculated based on the system 
operating conditions. Degradation mechanisms caused by physical, 
chemical or electrochemical reactions. These mechanisms affect 
the components characteristics such as the length of triple phase 
boundaries (TPB) and the electrolyte conductivity. Based on the 
values of degraded parameters, the system input/output (I/O) per- 
formance model calculates system performance indexes (such as 
system efficiency or output power). At the next time step, updated 
values are applied to the degradation model and this cycle is 
repeated until the end of the system operation lifetime. Up to this 
stage, the system performance deterioration as a function of oper- 
ating conditions can be presented. At the next step, using this func- 
tion, the optimum decision variables including operating 
conditions will be derived. 

Generally, optimization models tries to find the operating con- 
ditions or/and design parameters which minimize the cost of elec- 
tricity (COE), maximize the efficiency or profit through system 
operating lifetime |7]. In this study, our developed DBO framework 
optimizes operating conditions to achieve a minimum cost of elec- 
tricity. The DBO framework addresses the influences of a given 
operating condition on the degradation rate of the system and 
includes degradation effects in the techno-economic optimization 
[5]. 


T. Parhizkar, R. Roshandel/ Energy Conversion and Management 133 (2017) 20-30 21 


Nomenclature 


area, m? tortuosity 

gas phase concentration, mol m~? diffusion volume of simple molecules, cm? 
natural gas cost, $m? fraction of the reaction heat that is generated at the an- 
mass diffusivity, m? s7! ode 

activation energy, kJ mol”! pre-exponential factor, Q7! m~? 

Faraday’s constant, C equiv! overlap angle, 15° 

enthalpy, J mol”! bulk oxygen coverage, dimensionless 
current density, A m~? exponential activity parameter, kJ mol~! 
working current, A perimeter 

constant coefficient 

rate constant for oxidation reaction, m? kg! s~! Subscript and super script 

Nickel particles growth rate initial condition 


length, m bipolar plate (interconnect) 
natural gas mass flow rate, m? s~! fuel channel 


molecular weight, kg mol! anode 
equivalent electron per mole of reactant, equiv mol”! electrolyte 
number of nickel particles, mol cathode 
pressure, Pa air channel 
dimensionless pressure activation 
site occupation probability critical 
universal gas constant, J mol”! K~! channel 
Nickel particles radius concentration 
rate of oxidation, mol kg! s7! effective 
standard entropy, J mol”! K~! fuel 


time maximum value 
temperature, K i minimum value 


triple phase boundary, m m~? open circuit voltage 


decision variable Ohmic 
voltage, V optimal 
state variable oxydation 
molar fraction of component in the mixture reference 
Ohmic resistance, Q x direction 
over potential, V y direction 
electrical conductivity, Q`! m! z direction 
porosity 


S Qs DXX 





i I 
I 
i I 
I/O Performance |= = = > | Degradation | ı 
model f model 
| | 
L 


Degradation based optimization 
model 





Optimum operating conditions 


Minimum long term COE 


Fig. 1. Proposed framework for considering degradation mechanisms in optimization model. 


1.2. Background respect, recent studies aim at improving the understanding of 
degradation in SOFC stacks to manage system degradation rate in 

The relatively high investment cost of SOFC stacks leads most order to develop cost-competitive SOFC stacks applications in long 
studies to provide an analysis of possible scenarios that make SOFC term operation. The degradation is the process of systems perfor- 
applications feasible from an economic point of view [8]. In this mance deterioration; results from irreversible and gradual growth 
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of damage that happens during systems’ life cycle [9]. Degradation 
models are data driven or developed based on physical, chemical 
and electrochemical principals (principal based models) [10]. 

Nakajo et al. [11] presented a data driven model to investigate 
the effect of different operating conditions on the life time opera- 
tion of a SOFC stack. The developed model considers main degrada- 
tion mechanisms includes the reduction of electrolyte ionic 
conductivity, metallic interconnect corrosion, anode nickel parti- 
cles coarsening and cathode chromium contamination. The results 
show that lower stack temperature extends the life time of the 
Stack. 

Liu [12] developed a principle based model to provide the dis- 
tribution of thermal stresses during thermal cycling of a SOFC. 
The main contribution of the proposed model is to predict the life 
time of a sample cell. The results show that the current density and 
cell temperature are the most critical parameters that affect the 
SOFC microstructure. 

A data driven model based on neural networks is developed by 
Marra |13]. The electrolyte degradation is considered as the main 
degradation mechanism, and voltage reduction is then presented 
to highlight the cell deterioration. 

In a paper written by Nakajo et al. [14], a data driven model is 
presented to evaluate the probability of failure (PF) of a SOFC. The 
electrolyte degradation, metallic interconnect corrosion, anode 
nickel particle growth and chromium contamination of the cath- 
ode are the main degradation mechanisms which are considered 
in their model. Results show that adjustments of operating condi- 
tions, such as flow rates, temperatures, system specific power and 
methane conversion fraction in the reformer. These parameters 
decrease the failure probability of anode and cathode by a factor 
of around 300 and 30, respectively. 

The overall resistance model of a SOFC is developed by Aguilar 
et al. [15]. The sulfur poisoning is considered as the main degrada- 
tion mechanism. Impedance measurements show that the overall 
cell resistance can be reduced by diminishing the H2S content of 
the entering fuel. 

As can be concluded, the data driven models are more com- 
monly applied to analyze the degradation mechanisms. However, 
in principle based models, the effect of operating conditions on 
the rate of degradation mechanisms can be well determined | 10]. 
Moreover, the principal based models can be integrated in opti- 
mization frameworks which is the main outcome of this study. 

Many studies on SOFC stacks optimization have been done in 
the last decade [6,16]. Mostly the objective of the optimization is 
maximizing system profit or minimizing system cost regarding 
design or operating parameters as the decision variables |6]. Based 
on the reviewed researches, the approximation of system degrada- 
tion is a well-developed field and a large stream of research are 
conducted with various purposes. Similarly, operation optimiza- 
tion frameworks regardless considering components degradation 
are also available. However, obviously, there are not much litera- 
tures that combining these two fields. 

In this work, the SOFC system operating conditions are opti- 
mized regarding degradation mechanisms based on our developed 
DBO framework. The objective of present study is to minimize the 
COE of the SOFC stack through its operating life time. The key out- 
comes resulting from the proposed DBO framework are the stack 
optimum operating voltage and fuel utilization which lead SOFC 
stack to operate at the minimum COE over its lifetime. 

This paper is organized as follows: in Section 1, the evolution of 
the SOFC degradation and optimization models are presented. In 
Section 2, the main degradation mechanisms are introduced and 
degradation model and I/O performance model of the SOFC are pre- 
sented. Section 3 presents the optimization framework and the 
solution approach. Results are presented in Section 4 and finally 


the study is concluded with an overall discussion on DBO frame- 
work in Section 5. 


2. Modeling approach 


The general operating principle of a SOFC involves the trans- 
portation of oxygen from air at the cathode side to the anode side 
via O°% ions. At the anode side hydrogen is oxidized to form water 
and release electrons as shown in the reaction (a). The released 
electrons flow through the external circuit to the cathode side in 
order to produce electricity and simultaneously react with the oxy- 
gen to produce O° ions as represented in the reaction (b). It 
Should be mentioned that hydrogen is generated by an internal 
reforming process which takes place at the anode side. 


H + 0* —H,0+2e (anode) (a) 


O,+4e- — 207 (cathode) (b) 


The studied SOFC stack is 150 kW with the stationary applica- 
tion. The developed I/O performance model consists of six control 
volumes: the solid bipolar plate interconnects, the fuel channels, 
the anode layer, the electrolyte layer, the cathode layer, and the 
oxidant channels. The properties and constant values used in the 
simulations are referred to the Wen et al. [6] research. 

The main degradation mechanisms of SOFC stacks and the pro- 
posed degradation model of studied SOFC is presented in the fol- 
lowing sections. 


2.1. Degradation mechanisms 


During operation, SOFC stacks degrade in many different ways. 
A number of degradation pathways have been identified that can 
be categorized based on the occurrence part in the stack as pre- 
sented in Fig. 2 [17]. 


2.1.1. Degradation mechanisms in the interconnections 


e Oxide layer Formation 

The high reactivity of the Zirconium oxide towards La and Sr, 
apparent at the interface between Lanthanum strontium mangan- 
ite (LSM) and yttrium stabilized zirconia (YSZ), is the basis of the 
formation of LajZr,0O7 or SrZrO3. Oxidized structure deteriorates 
the catalytic activity of the cathode. The key parameter which con- 
trols the tendency towards the formation of foreign phases is the 
composition of LSM, and in particular, the amount of Sr. Manganite 
with a manganese oxide excess is proven to be more stable which 
is used in the most SOFCs | 18]. 


2.1.2. Degradation mechanisms in the cathode 


e Structure changes 
Operating at high temperature range of 900-1000°C may 
change the microstructural properties. This destruction diminishes 
triple phase boundary (TPB) and hinders reactants diffusion |19]. 


e Reaction with chromium 


In order to decrease the degradation of interconnecting plates at 
high operating temperature, plates are coated by a certain quantity 
of chromium. The input hot air flows on the surface of the metal 
and causes the transfer of volatile chromium oxide and hydroxide 
from the surface to the cathode. This mechanism will decrease the 
TPB length at the cathode side mostly in the cathode supported 
SOFCs [19]. 
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Fig. 2. SOFC degradation mechanisms. 


2.1.3. Degradation mechanisms in the electrolyte 


e Phase changes, impurities and dopant diffusion 
During the system operation at a hazard condition, phase 
change and impurities diffusion will happen in the electrolyte layer 
which can decrease the ionic conductivity | 18]. 


e Crack formation 


System leakage or mishandled on/off cycles cause the penetra- 
tion of the undesirable oxygen to the anode side. The existence of 
the oxygen at the anode side can lead to oxidation of the nickel 
particles. As there is a considerable difference in the molar volume 
of Ni and its oxide, this degradation mechanism leads to the forma- 
tion of internal stress. The internal stress causes tension in the 
electrolyte layer and crack formation. The permeation of gasses 
through the crack decreases system performance and accelerates 
other degradation mechanisms | 18]. 


2.1.4. Degradation mechanisms in the anode 


e Coarsening and oxidation of nickel particles 

The nickel phase in the yttrium stabilized zirconia (YSZ) struc- 
ture will be reformed under unsuitable operating conditions. 
Nickel particle size increases in this undesired reformation mech- 
anism. Moreover, nickel will be oxidized in the presence of oxygen, 
water, and oxygen ions at the anode side. The increase in the nickel 
particle size and nickel oxidation, diminishes the contact within 
the Ni phase. Therefore, electronic conductivity and triple phase 
boundary will be diminished [20]. 


e Crack formation 


As mentioned before, the oxidation of nickel can cause crack 
formation in the anode layer [18]. 


e Carburization on nickel particles 


During the internal reforming, hydrocarbons react with steam 
to form a synthesis gas containing carbon monoxide (CO) and 
hydrogen (H2). Then, the steam reacts with the carbon monoxide 
of the synthesis gas to form carbon dioxide (CO2) and more hydro- 
gen. This blend is rich in both (CO2) and (H2). In this mechanism 
carbon is deposited on the catalyst layer and as a result reforming 
process and electrochemical activity will decrease [21]. 


e Formation of nickel sulfide 


The presence of sulfur in hydrocarbon fuels leads to the forma- 
tion of Ni sulfide. This mechanism decreases the activity in the 
anodic catalyst |15]. 

The SOFC performance deterioration caused by all of these dif- 
ferent contributions. However, as can be concluded, in actual oper- 
ating conditions with pure fuel, the main degradation mechanism 
in anode supported SOFC stacks with Ni-YSZ and monolayer LSM 
cathode is nickel coarsening and oxidation which is considered in 
the present paper. 


2.2. Governing equations 


In this section, a principle based degradation model is devel- 
oped. As system works, nickel coarsening and oxidation occur 
which lead TPB and ionic conductivity deterioration. Different 
researchers used image treatment and analysis to study nickel 
coarsening mechanism. However, particle coarsening can follow 
the charging capacitor law which is suggested by Tanasini et al. 
[22]. 


rey = (r™™ — r°) [1 — exp(—Kgcapt)| + r° (1) 


where r° and r™* are the initial and maximum radius of particles 
and kscap is the rate of growth. In many studies the rate of growth 
is constant and it is derived by fitting to the experimental data at 
a desired operating conditions. In our study, the rate of growth is 
a function of temperature that is derived by fitting the experimental 
data in literature [22]. 

The second degradation mechanism which is considered in this 
study is nickel oxidation at anode layer. The main related reactions 
are presented in Fig. 3. 

As can be seen, nickel oxidation can occur via two different 
pathways. The first pathway is based on pure chemical reactions 
[23]. In this pathway, nickel can be oxidized in the presence of oxy- 
gen or water vapor at anode side as reactions (c) and (d). 


Ni + 50 a NiO (c) 
Ni + H20 = NiO +H; (d) 


The second reaction pathway is the electrochemical oxidation 
of nickel which is taking place at the surface between ion 
conductor and metallic nickel [23]. The oxidation reaction is as 
reaction (e). 


24 
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Fig. 3. Main degradation mechanisms of Ni oxidation. 


Ni +0% — NiO + 2e- (e) 


In all reactions, the nickel oxidation rate can be calculated as 
Eq. (2). 


Tox = Koxc(1 — 00)? 


(2) 


kox is the constant rate of reactions which is calculated as Eq. (3). 


kox = K exp (- a) exp ( is 00) 


RT 
Moreover, c is the mole fraction of O2, H2O and O* at anode 
side that should be defined. The amount of H2O and O” depend 
on input fuel and oxidant flowrates. However, oxygen can exist 
at the anode side because of leakage or water electrolysis as reac- 
tion (f) [24,25]. 


(3) 


(f) 


The oxygen partial pressure of this reaction can be calculated 
from Eq. (4). 


log (Poo), 150 =2 log Ga 


H2 


1 
H20 + H2 +502 


_ 26,000 


T +5.94 


(4) 

Nickel coarsening and oxidation reduce electrical conductivity 
and triple phase boundary [26]. The TPB length is calculated as 
Eq. (5). 


TPB(t) = 2[(m — 20)r(t)|twi(T) (5) 


r(t) is the nickel radius which is degraded due to nickel coarsening. 
In addition, following equation is used to evaluate the perimeter of 
nickel, Tyni, as a function of the number of sites in a cluster and the 
site occupation probability [26]. 
... {1 —pr(t) 50.5 

i(t) = + 3.6ni(t)~ 
Tni (C) nie) ( pr(t) ) T ( ) 
pr is the site occupation probability which is calculated as Eq. (7) 
[26]. 


(6) 


Ce (t) 


pr(t) = exp mgg] + 





(7) 


Furthermore, the degraded surface also reduces electrical con- 
ductivity. The electrical conductivity as a function of site occupa- 
tion probability is as Eq. (8) [26]. 


Geg (t) = Go(tan 6)[pr(t) — pr.]"~ (8) 


As discussed, TPB and electrical conductivity are two parame- 
ters which are degraded through time. TPB affects activation and 
concentration over potentials at anode side. In addition, electrical 
conductivity increases ohmic over potential at anode side. The 
effects of TPB and electrical conductivity reduction on over poten- 
tials can be calculated by applying degraded values of TPB and 
electrical conductivity to the SOFC I/O performance model. The 
governing equations of system I/O performance model which leads 
to power generation calculation are presented in Table 1, [27]. 


3. Optimization procedure 


The optimization objective is to minimize the cost of electricity 
(COE) through system lifetime regarding degradation mechanisms. 
The objective function is presented as follow, whereas satisfying 
constraints related to the system. 


Max]|u(t),x(t)] = Investment cost + Fuel cost + Maintenance cost 
l 7 Power generation 
(27) 


The investment cost is the cost of SOFC stack and fuel cost is the 
cost of consumed natural gas during given operating lifetime as Eq. 
(28). 

HT 
Fuel cost = J (Mngt) x Cng) -dt (28) 
t 

The maintenance/replacement cost represents the replacement 
cost of a stack cell in a given operating conditions. In most studies, 
total maintenance/replacement cost is considered as a constant 
value. However, in this study, system maintenance cost is corre- 
sponding to the rate of a degradation mechanism. In this way, sys- 
tem actual operating cost can be calculated continuously. The 
maintenance cost is calculated as Eq. (29). 


dni(t) 
dt 





t+T 
Maintenance cost = J ( x replacement cost -dt (29) 
t 
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Table 1 
SOFC I/O performance model governing equations. 


Output voltage 
Actual operating voltage 


Open circuit voltage 


—(Hy,0—-Hx, —$Ho, )+Tx a -450, ) —RT In | 


V(t) = Vow — (Ejnorm;(t) =F Neonc.3(t) + Nconc,5 + Nact3 (t) T fas hi j = 1,3,4,5 (9) 





0.5 


oy (10) 


Puy (Po, 





Voc = 
Ohmic over potential 
Ohmic over potential in regions 1, 4 
Electrical resistances of interconnect and electrolyte 


2F 





p4 = Lyla 
Ohmic over potential in regions 3, 5 





Area of active TPB regions 
Average diameter of electrodes parallel pores 


Nonmj = 1B; j=1,4 (11) 
L L L 

herea ied Gn (12) 
z (13) 

I(Aj; opt (t)—Ly Lz)Lj opt (t) I(Lj—Lj opt (t)) —_ 

Nonmj(t) = TA. ODL) F GLUT SIO’ j=3,5 (14) 

Aj opt(t) = (1 + 2L; ope(t)/Dpj)Lyz (15) 

Dpj = Dwj9;/(1 — 9;) (16) 


Concentration over potential 
Concentration over potentials at cathode 


Concentration over potentials at anode 


Water partial pressure 
Hydrogen partial pressure 


Oxygen partial pressure 


Diffusivities in anode and cathode 


Po, 2) 
Pu, Pu, o,TPe (t) 
Pu, o,TPs (t)Ph, 0 


— RT 
Neonc,5 = arn ( 


Nconc,3 (t) = ni In 


RT% hl 





RTW3 Li 





Pu, rep (t) = pfXn, 


Po, BP = Pair — (Pair — X02 Pair) €XP (at Balle 


( 
( 
Px, 0,7B(C) = PpXv,0 + O bL ( 
( 
( 
1.43x10 7T! (My, +My, 0)°° ( 


2FD3(t)03 LyLz 
RTWsLsI ) 





D3 = 


1 1 
Dp(2My, Myo)” (i, +3 0) 


1.43x10-7T"75 (Mo, +Mn, )°° (23) 


D; = 


Activation over potential 


2 


1 1y 
Pair (2Mo, Mn, yer (18, +VN, ) 








Activation over potentials Nace s(t) = RT sinh” | (; iz) j=3,5 (24) 
` a loj ’ ’ 
iti ; P P E; ; 
Exchange current densities io;(t) = YaTPB(t) (2) (72) exp (- ar) j=3,5 (25) 
Output power density 
SOFC output power density power(t) = V(t) xi (26) 


ono is the rate of nickel degradation and replacement cost is the 


cost of a cell replacement because of nickel degradation.Power gen- 
eration term in the objective function is the generated power 
through system operating lifetime which is presented as Eq. (30) 





t+T 
Power generation = i Power (t) - dt (30) 
t 


Power(t) = V(t) x i x A(t) (31) 


V(t) and A(t) are actual operating potential and active surface area 
over time. The optimization procedure is also subject to several con- 
straints. It is assumed that SOFC is designed to work at a specific 
temperature and current density ranges. The proposed constraints 
are expressed in Eqs. (32) and (33). 


Tmin < T < Tmax (32) 


jmin Z i 2 jmax (33) 


Moreover, it is assumed that the SOFC stack works until half of 
nickel particles in electrode is degraded [30]. In this case SOFC 
stack fails and cell should be replaced. 


0 
, ni 
Ae 
ni(t) 5 


Finally, the initial condition of oxidized surface area is equal to 
zero as Eq. (35). 


09 = 0 (35) 


(34) 


To solve the optimization problem, evolutionary techniques can 
be used to find the optimum decision variables due to their versa- 
tility and ability to optimize complex systems | 30]. In particular, 
the proposed optimization model has many local optima solutions 


with nonlinear objective function. Moreover, its parameters 
changes over time. Genetic algorithm (GA) technique appears to 
be a promising algorithm for handling this kind of optimization 
problems. 

GA algorithm has a parallel-population based search with 
stochastic selection of many individual solutions, stochastic cross- 
over and mutation that can lead model to the optimal solution in a 
short time. The GA algorithm initializes a population with random 
candidate solutions. Each candidate solution is called a ‘population’ 
which determines a vector of operating conditions. GA algorithm 
evaluate the objective function of the population and take its infor- 
mation into account to generate a new set of operating conditions. 
The population moves through the solution space following some 
basic formulae in search of a global optimum. After several gener- 
ations, only the ‘most optimum’ one can transmit information onto 
others, making the optimization very fast in comparison to other 
evolutionary techniques [31]. 


4. Results and discussion 


To evaluate the performance, effectiveness and efficiency of the 
DBO framework of SOFC stacks, it is applied to a 150 kW Ni-YSZ 
anode supported SOFC with a stationary application. The physical 
properties utilized to obtain simulation results are referred to 
Wen et al. [6] research. 


4.1. Model validations 


To validate the I/O performance model, a comparison between 
model outputs and results published by Wen et al. [6] is 
undertaken. The comparison illustrates that under the same SOFC 
geometric dimensions and operating conditions, the modeling 
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output accuracy is compelling enough to validate the numerical 
simulation. Furthermore, the voltage degradation rate at different 
current densities and fuel utilizations are used to validate the com- 
ponent degradation model. The results show that the model pre- 
diction agrees with the experimental data gathered by Muller 
et al. [32] and the results are nearly close. 


4.2. System characteristic 


Fig. 4 presents the V-I curves, at start and 10,000 h of system 
operation at three different temperature ranges. As it appears, 
the temperature affects the shape of the power profile. By increas- 
ing the temperature, the maximum power value moves to the right 
of the graph indicating the higher current density. Operation at 
high temperature range also results in low activation loss, which 
leads system to operate at a higher level of power and voltage. 
On the other hand, the degradation rate is found to increase with 
higher temperature. It is worth noticing that high temperature 
ranges, which correspond to higher voltage and output power, lead 
to an acceleration of the degradation and result in higher voltage 
and output power deterioration rates. In other words, in order to 
reach higher output power in long term operation, power deterio- 
ration rate dependency of temperature should be taken into 
consideration. 

Fig. 4(c) shows that the voltage is degraded about 12% after 
10,000 h of operating at a high temperature range (1108 K). The 
breakdown of SOFC overpotentials after 10,000 h of operation is 
presented in Fig. 5. 

As can be seen in the figure, the activation overpotential which 
originates from irreversibilities of electrochemical reactions has 
the greatest portion of overpotentials. 


1.4 250 
S ir 200 = 
© 08 150 = 
FA 0.6 100 $ 
o 0.4 6 
> 02 90 a 

0 0 
5 





Voltage (V) 





(b) 


1.4 250 
12 =~ 

200 
=. = 
o 0:8 al 
g o6 100 Ẹ 
= 04 en, S 
> 0.2 a 

0 0 





Current density (A/cm?) 


(c) 


Fig. 4. Polarization and power curves at start (—) and after 10,000 h (---) for (a) 
low-1048 K (b) medium-1078 K (c) high-1108 K temperatures. 


In contrast, the concentration overpotential that is generated by 
mass transport in electrodes, is often neglected due to the high gas 
diffusion rate at high operating temperatures. Therefore, as is 
shown, it has the least value [33]. The last overpotential is the 
ohmic over potential which is generated by the ionic and electronic 
charge transfer resistances. The ohmic resistance is important in all 
types of cells and is essentially linearly proportional to the current 
density. For an anode supported SOFC, the ohmic loss is smaller 
than in the other types of SOFC cells [33]. Results indicate that 
around 28% of the total over potential is due to the ohmic loss. 

Each overpotential is caused by system inevitable physical, 
chemical and electrochemical losses (actual losses) and losses 
which are added to the overpotentials because of system degrada- 
tion mechanisms (degradation losses). 

Fig. 5 depicts the proportion of these losses on each overpoten- 
tial. For instance, a quarter of the activation loss occurs as a result 
of degradation mechanisms. The TPB is the most important param- 
eter which affects the activation overpotential. The TPB length is 
mainly affected by the microstructure of the electrodes. Nickel 
coarsening and oxidation decreases TPB length drastically. 
Additionally, based on the results, degradation has the most effect 
on the concentration overpotential, and about 45% of the concen- 
tration overpotential is due to the anode degradation. That is 
because of the fact that nickel degradation mechanism changes 
the structure of materials and decreases diffusivity dramatically. 
Furthermore, at a high temperature level the electrical conductiv- 
ity of nickel oxide and nickel particles are nearly the same. 
Consequently, the effect of nickel degradation on ohmic loss is 
negligible. 


4.3. System degradation 


In this section, the outputs of the degradation model are 
reported. All the results are obtained under the similar geometrical 
features and the same operating conditions reported in Section 4.2 
(47% fuel utilization at three different temperature ranges). Fig. 6 
(a) depicts the behavior of the active surface area over time. While 
the system is working, nickel coarsening and oxidation occur. The 
increase in the nickel particle size and nickel oxidation, decrease 
the contact within the Ni phase and, therefore, the active surface 
area is diminished. As nickel coarsening and oxidation rates are 
accelerated at upper temperature ranges, the deterioration of 
active surface area would be higher. 

Fig. 6(b) indicates output power over time. This figure shows 
that the output power decreases through time as a consequence 
of TPB and anode layer conductivity reduction. Furthermore, after 
a while, the power value becomes constant, depending on the 
operating temperature. In addition, as can be seen, at upper tem- 
perature ranges, the initial output power of SOFC is higher, but 
the deterioration rate is accelerated as well and power falls drasti- 
cally through time. On the contrary, at lower operating tempera- 
tures, the system operates at a lower power level with a smooth 
deterioration rate over time. 


4.4. System productivity and cost of electricity (COE) 


As discussed in Section 4.3, operating conditions have a great 
impact on the system power generation and efficiency in a long 
term operation. In this section a voltage and fuel utilization sensi- 
tivity analysis for system productivity and cost of electricity is per- 
formed. In Table 2, the economic assumptions applied to the COE 
calculations are summarized. 

The influence of the cell voltage on power generation through 
system lifetime and cost of electricity is shown in Fig. 7(a). The 
productivity reaches its maximum value at 0.51 V and the system 
generates 5856 kWh power through its operating lifetime. It 
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Fig. 5. Degradation based and actual losses after 10,000 h operating at (1108 K, 4 A/cm?). 


should be noted that the COE shows its minimum value at the 
same voltage. At low cell voltages, output power is low, and it 
causes power generation to be at its minimum level. As cell voltage 
increases, system operates at upper power levels and the power 
generation increases till it reaches its peak. By increasing the cell 
voltage, the degradation rate is accelerated, therefore, system oper- 
ating lifetime decreases and system performance deteriorates 
more rapidly. This lead system to have a lower productivity (power 
generation over the system lifetime). 

Moreover, the sensitivity of COE is illustrated in Fig. 7(a) as well. 
At low cell voltage the electrical plant efficiency is low. As a conse- 
quence, a large natural gas stream is consumed, which is a cost 
dominant factor. Moreover, as mentioned, at low voltage, the pro- 
ductivity is low. High fuel cost and low productivity increases COE 
at low voltages. At higher voltage ranges, the system becomes 
more efficient and generates more power until the COE reaches 
its minimum value. At higher cell voltages, system overall effi- 
ciency deteriorates drastically, therefore, fuel cost increases and, 
as before mentioned, system productivity decreases. These facts 
cause an upward trend in the COE. 

In Fig. 7(b), the sensitivity analysis of fuel utilization is illus- 
trated. Fuel utilization is related to the natural gas feed stream. 
As fuel utilization increases, system power density rises and, as a 
result, system productivity increases up to 3000 kWh, provided 
when fuel utilization is adjusted to about 60%. At higher fuel uti- 
lization ranges, system operates at upper power densities but the 
degradation mechanism rate increases and performance deterio- 
rates more rapidly as well, thus productivity drops. 

Fuel utilization variation has the same effect on COE trend as 
voltage variation. Increasing fuel utilization, enhances system effi- 
ciency and productivity until system reaches its minimum COE 
value. At higher fuel utilization, performance and efficiency deteri- 
orate drastically because of the higher degradation rates and this 
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leads system to have a higher fuel cost with lower productivity. 
These two effects cause system to operate at higher COE. 


4.5. System optimization 


As already discussed, SOFC operating conditions have a great 
impact on system productivity and COE. The optimum operating 
conditions which lead system to have a minimum COE are deter- 
mined according to the proposed degradation based optimization 
(DBO) model. To study the advantages of the DBO model, two opti- 
mization strategies are compared. In the first strategy, system 
operating conditions are optimized to maximize power density 
regardless of degradation mechanisms. Usually, these operating 
conditions are defined by the manufacturer of the system and here 
in it is based on Wen et al. [6] optimization results (base case). In 
the second strategy, system operates at the conditions which are 
derived from our proposed DBO model. The derived optimum 
operating conditions are presented in Table 3. 

Fig. 8(a) presents the predicted COE of the system at mentioned 
operating strategies (Table 3). As can be seen, operating under the 
DBO strategy, decreases system COE up to 7.12%. Under the base 
case strategy, the system operates at a higher power level, but 
degradation mechanisms are not controlled and performance dete- 
riorates drastically. Moreover, uncontrolled degradation mecha- 
nisms, cause system to fail earlier, therefore, the power 
generation decreases. On the other hand, the DBO strategy, consid- 
ers degradation mechanisms. This leads system to operate at a 
lower power level over a longer operational lifetime, therefore, sys- 
tem COE would be lower. The breakdown of the system lifetime 
COE for both operating strategies is illustrated in Fig. 8(b). This fig- 
ure indicates that the proposed DBO strategy leads to a more effi- 
cient and durable operation of SOFC stacks. At the DBO strategy, 
the power deterioration rate is much lower than in the base case 
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Fig. 6. Cell active surface area and output power over time at three temperatures. 
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Table 2 
Assumed economic parameters. 


Parameter Value Reference year 
Investment cost 1620 $ m~? [28] 

Interest rate 10% - 

Fuel cost 0.03 $ kW h™! [28] 

SOFC Single cell cost 90 $ kw! [29] 


and system operating lifetime is longer. As a result, system gener- 
ates more power. It should be noted that most of the SOFC lifetime 
cost is in the initial investment of the system, which is the same for 
both operating strategies. Furthermore, although power generation 
is lower in the base case, the fuel cost is higher in comparison with 
the DBO strategy. The main reason is that the long term efficiency 
in the base case is much lower than in the DBO due to the acceler- 
ating degradation rates. 


4.6. The economic parameters’ effect on the optimum COE 


The contributors to the system cost are the investment of the 
SOFC stack, the maintenance/replacement and the fuel cost. The 
notable fraction of the total cost is in the investment. Conse- 
quently, the primary focus of the total cost reduction should be 
centered on SOFC stacks. In recent years, new materials and high 
manufacturing technology have sharply reduced the SOFC price 
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Fig. 7. Voltage and fuel utilization sensitivity analysis of COE and productivity. 


Table 3 
Optimum operating conditions resulting from DBO. 
Value Unit 
Fuel utilization 40 % 
Voltage 0.60 V 
Cell temperature 1033 K 


and The Fuel Cells and Hydrogen Joint Undertaking (FCH JU) has 
predicted that the stack cost is projected to be 2000 $/kW by 
2020 [34]. In addition, the fuel cost variation can have different 
effects on the COE. The sensitivity of the COE to the reduction of 
investment cost to 2000 $/kW and the fuel cost fluctuation of 
+20% is performed in this section. The results illustrate that system 
economics could be improved by decreasing fuel and investment 
cost. However, COE is less sensitive to fuel cost compared to invest- 
ment cost because the main portion of the system cost is in the ini- 
tial investment. 

At the current fuel cost, decreasing system capital cost by half 
could potentially reduce the COE by around 30%, and the SOFC sys- 
tem cost of electricity can be more competitive with the grid. In 
addition, when the investment cost is lower, the impact of the fuel 
cost on COE will be greater. For instance, 4% COE reduction is 
obtained by a 20% reduction in fuel cost at the current investment 
cost. However, by cutting the investment cost in half and reducing 
fuel cost by about 20%, the COE could decrease up to 6%. 
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Fig. 8. Comparison of COE, productivity and fuel cost for optimum and base case 
conditions. 


5. Conclusion 


This research aims to contribute to the field of degradation 
based optimization methodology. The proposed model optimizes 
operating conditions of an anode supported SOFC system regarding 
degradation mechanisms to achieve minimum cost of electricity in 
long term operation. Based on the suggested degradation model, 
the influence of different operating conditions on anode side 
degradation and its effect on performance deterioration is illus- 
trated. Moreover, the effect of operating conditions on power gen- 
eration and cost of electricity is also determined and the optimum 
operating conditions to achieve minimum COE are derived. To 
show the effectiveness of the proposed framework, results are 
compared with the optimum base case strategy which its operat- 
ing conditions are optimized regardless of degradation mecha- 
nisms. Results show that using the DBO framework, system COE 
is 7.12% lower in comparison with the base case. Moreover, the 
sensitivity analysis of COE indicates that the COE is less affected 
by fuel cost than investment cost. 

In this study, the validation and demonstration of the degrada- 
tion based optimization model laid the framework for considering 
degradation in long term performance analysis of SOFCs. However, 
this framework can be expanded to any other energy systems 
including power plants, and refineries. 
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